Spatial correlations of singularity strengths in multifractal branching processes 
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The n-point statistics of singularity strength variables for multiplicative branching processes is 
calculated from an analytic expression of the corresponding multivariate generating function. The 
key ingredient is a branching generating function which can be understood as a natural generalisation 
of the multifractal mass exponents. Various random multiplicative cascade processes pertaining to 
fully developed turbulence are discussed. 
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I. INTRODUCTION 



Many complex processes occurring in nature are believed to be organized in a selfsimilar way: turbulence, large- 
scale structure formation in the universe, high- energetic multiparticle dynamics and diffusion limited aggregation are 
but a few examples showing scaling over a wide range. 

A convenient tool to characterize the spatial configurations of these processes is the multifractal formalism [Q- H 
which analyses the one-point statistics of singularity strengths under the assumption of multiple local scale invariance. 
Typically, scale-invariance is a consequence of some hierarchical organisation of the underlying process. Since this 
process occurs in space, its characterisation in terms of multifractals can be expected also to yield information about 
. spatial correlation functions. 

Multivariate correlations, on the other hand, have long been known to contain (at least in principle) complete 
information of the underlying random process. Moments and especially cumulants are capable, in their multivariate 
form, of providing this information by means of well-tested statistical procedures. As in the case of multifractals, 
i-j^"* cumulants can detect and characterize scale invariance, but unlike the former do not depend on its presence. 

The existence of these two sets of tools naturally leads to the question how they are interrelated: given a selfsimilar 
process, does the multifractal characterisation completely determine the spatial correlation functions and vice versa, 
or does one provide more information than the other? 

Some attempts have been made to elucidate the relationship between multifractals and correlation functions: two- 
point statistics of multifractal measures have been discussed in Refs. [^,|| and applied to turbulence in ||. It is to be 
expected, though, that only the full n-point statistics of singularity strengths would provide the equivalent information 
contained in the spatial correlation functions of all orders. 

It is known that different selfsimilar processes can lead to identical multifractal exponents [0-|9|. This can be 
understood as a hint that the multifractal characterisation is indeed incomplete and that the spatial correlation 
functions may contain more information about the underlying process. We shall show that this is indeed the case for 
several examples of multiplicative branching processes. 

The ultimate goal for any analysis, of course, would be a characterisation of the multivariate statistics gained 
from a physical process to a degree of completeness that would permit precise reconstruction of the "branching rule" 
that governs each step of the selfsimilar process. As far as theory is concerned, this holy grail would be achieved 
with the specification of a multivariate generating function whose n-fold derivatives would precisely reproduce the 
corresponding experimentally measured n-variate cumulants. 

Clearly, finding this real-life generating function remains beyond present-day capabilities, for both theoretical and 
experimental reasons. It is possible, though, to make progress towards the ideal by inventing simple branching models 
whose generating functions can be calculated. An approach to calculate the n-point spatial correlations to arbitrary 
order within some of these models was presented in Refs. [[lO 11 1: there, the multivariate generating function of the 



spatial correlations was constructed iteratively from a backward evolution equation, leading to a recursive derivation 
of spatial correlations. An important further step has been the recent discovery of a large class of analytic solutions 
for generating functions of selfsimilar multiplicative branching processes fi~2|| . 
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In this paper, we expand on the latter discovery and attempt to cast further light on said relationship between 
multifractals and correlations. After a brief explanation of the mechanics of multiplicative branching processes in 
Sec. H we review in Sec. pj the formalism of multivariate correlations and derive the abovementioned analytic 
solutions for the branching generating functions. While singularity strength variables are already treated in Sec. Ill, 



the relationship between spatial correlations and multifractals is explored more fully in Sec. IV, where we also show 
how the vaunted splitting function ("branching rule") can be reconstructed via a two-dimensional Laplace transform. 
The question of reconstructing this splitting function from data is also briefly discussed. Four examples of branching 
processes, treated in Sec. [v| drive home the message that multifractals contain less information than the spatial 



correlations. We discuss our results in more general terms in Sec. VI 



II. MULTIPLICATIVE BRANCHING PROCESSES 



Random multiplicative branching processes, sometimes also called weight-curdling models or simply cascade models, 
can be constructed for any number of branches per splitting. For simplicity, we concentrate on binary processes; 
generalizations are straightforward. 

The branching rule governing a binary cascade is described as follows: a starting energy density e^ ^ = 1, uniformly 
distributed on the unit interval, is split up into two daughter densities — go e ^ and = qie^ living on adjacent 
subintcrvals of length 2~ 1 . The random weights q$ and q±, often also called multipliers or splitting parameters, are 
drawn from a joint probability density (or "splitting function") p(qo,qi). In the next generation, each of the two 
daughter energy densities is itself split up by the same branching rule into two granddaughters distributed uniformly 
over adjacent subintervals of length 2 -2 . Generally, energy densities of the j-th generation, eW, are characterized by 
the binary index k — (fcife ■ • ■ kj), with each k taking on possible values or 1. Successive repetition of this prescription 
for each e of the j-th generation yields, at the next branching, 2 J+1 energy densities ^^k-o = ^11— fe-o e fe?---fe- anc ^ 

e k*-k-l = Ik^k-^ki—k-i with the multipliers drawn from a splitting function K'Zfe^ifcO' Qk^ ^ l)' which usually 
(although not necessarily) is identical with p(qo,qi). The process is completed after J cascade steps. 

When the splitting function p(qo, qi) does not depend on the generation j or branch location k\ ■ ■ ■ kj, the above 
prescription gives rise to a selfsimilar multiplicative branching process, wholly characterized by the choice of splitting 
function p. We list four examples for later use. The choice 

p(q , gi ) = I [% - (1 + 0)) + S(q Q - (1 - 13))] S(q a + q x - 2) (1) 

leads to the binomial multiplicative process, also known as the p-model (often the equivalent parameter p = (1 + 0)/2 
is used) 0. The delta function 5(qo + q\ — 2) ensures that this process conserves energy at every branching. 
The splitting function for the non-energy conserving counterpart of the p-model, the a-model ||, is given by 

l 



p(«d,?i) = ~II + /?)) + % fc -(l -13))] . (2) 



fe=0 

For the energy dissipation process in fully developed turbulence, the parametrisation 

lT(2/3) (q Q \P-^ (q x \P 



Kfc.fc) = 2fo8)ah-j It) 2) 



1 Tm (q Q (2-q ) f- 1 6(q + q 1 -2) (3) 



2 2 / 3 - 1 r(/3) 2 



with (3 = 3.2 has been directly deduced from experiments p3[ . The corresponding splitting function not conserving 
energy, is given by 

P(go, Qi) = ( 22j L r^ ) 2 (go(2 - q^f- 1 ( gl (2 - q^ 1 . (4) 
These four splitting functions are shown in Fig. 1. 
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FIG. 1. Splitting functions p(qo, qi) of (a) the p-model with = 0.4, (b) the a-model with = 0.4, (c) the energy-conserving 
SRST cascade model with = 3.2, and (d) the energy-non-conserving SRST cascade model with = 3.2. 



III. SPATIAL CORRELATIONS 



A. Possible variables 

For branching processes, the success in extracting useful information or presenting a clear picture depends strongly 
on the choice of variable. In Refs. Jl^,|ll|,0, forward and backward evolution equations for a multivariate generating 
function were employed to calculate the correlations directly in terms of the energy densities, {tj^l • ■ In the 

context of multifractals on the other hand, another set of variables, the so-called singularity strengths aj^ , are 
in use. The definition of these strengths, (2~ J )e^ = (2~ J ) a K , means that they are related to the energy density 
variables e$ by 



„( J ) 



(J) 



In 

Jln2 n 



(5) 



so that correlations (a^ ■ ■ ■ ) between the singularity strengths are intimately related to corresponding 



correlations ((lne^) • • • (lne^)) between the logarithms of the energy densities e^ J . 



(J) 
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We have shown previously |12j that correlations between the lne lead to a particularly transparent structure for 
multiplicative branching processes. Starting with correlations of the logarithms of energy densities, we show below 
how this transparent structure arises before translating these results into the language of multifractals. 



B. Correlations and generating functions 

We consider the spatial correlations 

pK lt ...,K n 



In ■ • • In e K 



(J) 



(6) 



of order n between the variables rne^p. The bracket (. . .) indicates the averaging over all possible configurations. 
Moments Pk 1 ,...,k„ are conveniently calculated from the multivariate generating function 



Z[X^]= exp £ A« fcj lne« fcj 



(7) 



by taking appropriate derivatives with respect to the conjugate variables A^. 

d n Z[\^] 



(J). 



PKi,...,K„ 



A=0 



The corresponding multivariate cumulant generating function 

iqA (J) ] = lnZ[A< J )] 
yields the cumulant correlation densities in the same way via 

d n K[\^] 



lne (J) ---lne (J) 



(8) 



(9) 



(10) 



A( J )=0 



here, we have introduced the index c (for "cumulant") on the last bracket to distinguish ( |l0| ) from the expression 
for the ordinary correlation densities. The lowest four orders of cumulant correlation densities arc 



C Kl = (lne Kl ) c = (hxe Kl ) , 
Ck, u k, 2 = (lie/Ci hie K3 ) c 

= (In e Kl In e K2 ) - (In e Kl } (In e K2 ) , 
= (lne Kl lne/t 2 lne K3 ) c 

= (In e Kl In e« 2 In e Ka ) - ^ (In e Kl In e K2 ) (In e Ks ) 

{3} 

+2(ln e Kl ) (In e K2 ) (In e Ka ) , 
c Kl = (lne Kl lne K2 lne K3 hie K4 ) c 

= (lne Kl lne K2 lne K3 lne K4 ) - ^(lne Kl lne Ka lne K3 )(lne K4 ) 

{4} 

- ^ (In e Kl In e K2 ) (In e K3 In e« 4 ) + 2 (In e/cj In e« 2 ) (In ck 3 ) (In e K4 

{3} {6} 

-6(lne Kl )(lneK 2 )(lne K3 )(lneK 4 ) , 



(11) 
(12) 

(13) 



(14) 



with brackets {. . .} indicating summation over permutations of indices. 
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C. Construction of the cumulant generating function (first approach) 



The derivation of the cumulant generating function depends crucially on a rearrangement of terms in the sum 
Y)f !l kj ^%^...kj me fcf...fc, entering the generating function in (^). The trick lies in recognising that the cascade 
prescription of Section O implies that the energy density in the J-th generation is the product of all the splitting 
parameters in its ancestry, 



(J) _ (J) _ (1) (2) (J) 
e K ~ e kx---k ., — Ik! 9feife 2 '" Sfci- -kj ' 



so that its logarithm is additive in the multipliers' logarithms, 



i=i 



Inserting this into the sum in the exponential of the generating function in (Q) and defining 



. = V a ( ; 7) 

we find, after rearrangement and judicious regrouping of terms, that 

^ ^ki--kj me fei--fej = E E ^k-L—kj m 5fei...fcj 



J=l ki,...,kj 
J 



(15) 



(16) 



(17) 



(18) 



This can now be utilized to rewrite the configuration average (0(e^ J ')) of any function O in terms of the splitting 
parameters as follows. We start with the definition 

(0(e^)) = |d(ln^ ) )...d(ln^ ) Jp(ln e « ,...,ln e ^ ) 1 )O( e ( J )) , (19) 

where p(ln £q^ , • • • , In e^.i) is the joint probability of finding a given logarithmic energy density In £q^ in bin (0 • • • 0), 
a given lne^^ in the next bin, etc. Since each energy density is completely determined by its J ancestral splitting 
parameters, the multivariate probability density p(lneQ^ , • • ■ , lne^^) can be constructed from the splitting function 



p(q ,qi) by using (|16|): 

Klne« Q ,...,lne« )= f[f[ f[ dg£>... fcj _ l0 dgff..^ pfog>..*,_ l0 , gg..^) 



J 1 



j = l *1, 

fc j-l=° 

1 



fcl,..., J=l 



(20) 



Inserting (l8f), (19) and ( |20| ) into Eq. (|9j), we find that the cumulant generating function becomes 



K[\i J ! ,...,XY.> 1 ]=ln 



(J) 



d(ln e £. fcj ))p(lne« ,...,lne« 1 ) 



j i 



E E o[Ag 

j=i fc i 

fc,-„, =0 



U) x (i) 1 

fc,-_i0> A fci---fe,_ilJ ' 



(21) 
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where 



Q[Ao, Ai 



In 



dq dqi p(q , 9i) exp (Aq lng + Ai lngi) 



ln(exp (A In g + Ai lngi) 



(22) 



Eqs. (21 )— (J22)) represent the long-sought analytic expression for multiplicative cascades: they show that the cumulant 
generating function of the entire cascade can be written as the sum of all branching generating functions Q, one for 
every branching. In most cases, the two-fold integrals ( p2| ) can be solved analytically or parametrically to yield a 
complete analytic solution for K and thereby for all cumulants. 

The scope and limitations of the solution ( [2l| ) are as follows. Clearly, the generating function in terms of In e is 
applicable to any functional form of the splitting function or b.g.f. It does not depend on the number of branches 
either: trivariate or even higher-variate splitting functions can be implemented. Due to the additive nature of the 
b.g.f.'s, the splitting functions can differ from generation to generation and even from branch to branch. The only 
(and important) precondition for the applicability of Eq. (|2l]) is that the splitting variables of every branching must 
be independent of those of the other branchings and generations. 



D. Construction of the cumulant generating function (second approach) 

The cumulant generating function can also be constructed from a forward evolution equation. In Ref. pM, the 

(i) — 

forward evolution equation for the generating function was derived in terms of energy densities e)^ terms 
of the lne^...feo the cumulant generating function after j cascade steps, K^[X^], can be expressed in terms of 
IfO'-i^O'-i)] as follows. At each of the 2 J_1 independent branchings, Eq. ( |l6| ) relates the daughter energy densities 
to their parent by 



it. k] (% = o,i). 



In analogy to (|T^) , 



,0) 



fci,. 



n exp 



0) 



+ A 



(/) 



fci,..., 

Is < _ 1 =0 



ki—kj-iO A fci---fcj_il 
0) ,„ „0) 



In e 



k-i ■■■k-i 



I,, „0) I \0) 1 — „0) 



„0) 



and assigning, as in ([170, 



A 



0) 



A fel'"-fcj-l0 + Ajf l '...j(5 i _ 1 l — A fe 1 ---fc J _ 1 ! 



0) 



A 



0-i) 



we arrive at 



E ln [/ ^kl-k^^kl-k^i p^kl-k^^kl-k^i) 

1=1,..., j 

i f \W In „0) i \0) i„ „0) 

' ^ A fci,,,fcj_io Ui yfei---fe J _i0 T A /ci---fc i _ii in yfei---fc J _ii 



x exp i 



^o-D [A o-i) ]+ 53 g>w..^ i0! A«.. fcj ; ,; . 



(23) 



(24) 
(25) 



(26) 
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This forward evolution equation can be iterated from j — J down to j — 0, where we arrive at K^ =0 ^ [X^] = 
Ib.(Z^[X^}) — 0; the outcome of this iteration is identical to the solution (|2l|). Similarly, a backward evolution 
equation for the cumulant generating functions can also be used to achieve identical results; we leave this as an 
exercise to the reader. 



E. Cumulant correlation densities 



Having found an analytic expression of the multivariate generating function iffA^], it is now straightforward 
to calculate the cumulant correlation densities (|ll])-(|l4|) between the In t^) -variables via Eq. (|l0|). In order to 
give compact expressions, we define an ultrametric distance between bins. Assume the two bins are together for 
the first j steps of the binary cascade before splitting at the (J + l)-th generation, k% — (fci • • ■ kjkj + \ ■ ■ ■ kj) and 
1^2 = (fci • • ■ kjk'j +1 ■ ■ ■ fcj), where hi = k[ for all 1 < i < j and fcj+i ^ kj+i- Then the ultrametric distance between 



them is 



e?2 = dist(«i, k,2 
For n bins, the generalised ultrametric distance is 



J-J 



max dist(/Ci, Kj) 

l<i<j<n 



(27) 



(28) 



In the following, we assume that the same functional form for the b.g.f. is used at all branchings. We also assume 
the splitting function p(qo,qi) to be symmetric under exchange of the splitting parameters, p{qo 7 qi) = p{qi,qo)', the 
four splitting functions given in Eqs. (0)-(Q) fall into this category. In consequence, the branching generating function 
Q[Ao, Ai] = Q[Ai, Aq] is also symmetric, so that also 



9" 1+ " 2 g[A ,A 1 



dx^ax? 



d^+™Q[X ,X x 



A=0 



dX^dX? 



(29) 



A=0 



This relation results in compact expressions for the cumulant correlation densities C*Ki, Inserting ( pi|) into Eq. 
( |To| ) and taking into account (|l7|), the cumulant correlation density of first order is found to be 



Ck x — J 



( <9Q[Ao,A! 



V 9A 



A=0 



For second order we get, with (p7|), 



d 2 Q[X 0l X 1 



dXl 



+ (1 - *efa,o) 



d 2 Q[X Q ,Xi 



a=o/ \ dA <9Ai 

while with the generalisation (Eq), the third order cumulant density becomes 



A=0 



r (t 4 \ ( ^Q\MM \, h x / d 3 Q[X ,X 1 

C Ku k 2 ,k 3 -(J- da) ^ dX 3 ) + (1 - *«fe,o) [ dX 2 dXi 



A=0 



(30) 



(31) 



(32) 



In the case of fourth order, two cases have to be distinguished, depending on whether the first splitting of the four 
bins K\, . . . , K4 goes into three and one, say (/C1K2W3IK4) with dafoi, K2, K3) < c?4 and n\ = 3,?i2 = 1, or into two 
and two, say {k\k-i\k-$k\) with di(K\^Ki) < d±, g?2(k3,^4) < d^ and ri\ =2,n2 = 2: 



CKl,K2,K3,K,4 ~ {J _ ^4) 



9 4 Q[A ,Ai 



dXi 



A=0 



+ (1 - <$d 4 ,o) 



d 4 Q[X Q ,X 1 



OX^dX^ 



A=0 



(33) 



Even under the stated assumptions, the expressions (J30])-(|33D are still very general for selfsimilar binary multiplicative 
cascade models. The one and only input is the splitting function p(qo,qi) determining the branching generating 
function Q[Ao,Ai]. The cumulant correlation densities for the splitting functions ([j]) _ (|]) will be discussed in more 
detail in Sec. ^. 

It is clear from Eqs. (pfj|)-(p3|) that the derivatives of the branching generating function completely fix the spatial 
cumulant correlation densities. Using the definition ( p2| ) these derivatives can be expressed in terms of generic 
branching moments: 
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In the lowest three orders, they read explicitly: 



= ((\nq r(\n qi D c 



(34) 



A=0 



( ln go) c 
((lng ) 2 ) c 
((\nq a )(\n qi )) c 

((lng ) 3 ) c 
((lnq ) 2 (lng 1 )) c 



(In go) , 

((ln (Zo ) 2 )-(ln (? o) 2 , 
((hig )(higi)) - (lng ) (lngi) , 

((In go) 3 ) - 3 ((In go) 2 ) (In go) + 2 (lng ) 3 , 

((]xiq f(ln qi )) - 2 ((lng )(lngi)) (In go) - ((lng ) 2 ) (higi) 

+2 (In g ) 2 (In gi) . 



(35) 



F. Link to singularity strength variables 

We have defined the cumulant correlation densities Cf^, ...,/«„ — (hie^ • • 'lne^ ) c in terms of the In e^p- variables. 
With the help of (|), they can be transformed into cumulant correlation densities C£ ... n n — ( a K? ' ' ' a K,l)c between 
the singularity strength variables a~^} . We find for the corresponding generating function 



kj = 

= E A fei---fej + 



■lne^ 



fei •••kj 



(-i) 

Jln2' 



X {J) 



(36) 



Taking derivatives, the relations between the cumulant densities are found to be, for first order, 



and, for n > 2, 



CI 



i 



(Ine^) 



Jln2 



Jln2 



Cm 



-1 



(lne«...lne^) c 



-1 



Jln2/ ^ ' w V^ 1 n2 / 

in other words, with the exception of first order, the C% lt are directly proportional to C^j, ...,«„ 



(37) 



(38) 



IV. SPATIAL CORRELATIONS VS. MULTIFRACTALS 
A. Multifractals 

One way to approach multifractals is to count how often specific values a of the singularity strengths occur; this 
leads to the /(a)-spectrum jlj. The latter is related by a Legendre transformation to a sequence of mass exponents 
t(v) which are defined by the moment scaling behavior of bin energies E^... k .\ 



E 



E 



U) 

ki 



E 



eg) k 

23 



(39) 
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In order to avoid confusion we have written t(u) instead of the more familiar notation r(g). Since the bin energies 

are related to the singularity strengths by = (2~ j )"k , these moments are to be understood as scale-dependent 
measures for the one-point statistics of the singularity strengths. 

We first derive the expression for the exponents t(u) for the special case of a symmetric splitting function 
p(<70i9i) = p(qitQo) of a binary multiplicative cascade process which conserves energy. (Note that for energy non- 
conserving splitting functions, as for example those given in (^) and (^), the multifractal exponents cannot be defined 
unambiguously; see Ref. fll|| .) Then we can use ( (To]) and the statistical independence of splitting parameters to deduce 



') = K)KJ ■'•(&...!,,) = (<£)* ■ (40) 



Insertion into Eq. yields: 



V ; La 2 V2"- 1 / V ' In 2 V ' 



For a qo/'Zi-asymmetric, energy-conserving splitting function we have to replace the sum appearing in (39) by 



which gives rise to 



r 



E (KU)") = h [«> + «>]' • < 42 > 



w ' ,„(M>±M>u_ y+ '° m+m . ,43, 



In 2 V 2 " J ln2 



B. Relationship between Q[Ao,Ai] and t(u) 

The close relationship between the branching generating function Q[Xo, Xi] and the multifractal exponents t(v) is 
revealed by simple inspection of (123) and (E3J) . We find that 



r(z/) = -i/+— ln[exp(Q[A = ^Ai = 0]) + exp(Q[Ao=0,Ai = i/])] . (44) 

This relationship^ simplifies further once we consider an energy-conserving, qo/qi symmetric splitting function. 
Through Q we arrive at 

t(i/) = -(i/-1) + ^Q[Ao = i/,Ai = 0] . (45) 

The two relations ( p[ ) and ( |45| ) express t(v) in terms of Q[Ao, Ai]. It appears that r is more limited than Q in that 
the latter is defined over the full Ao-Ai plane while t(v) is restricted to the axes Ao = or Ai = 0. 

Let us consider the point in more detail. In particular, we will demonstrate that different energy-conserving splitting 
functions may lead to identical multifractal exponents while their branching generating functions do differ. 

Given an energy-conserving, go ll\ -symmetric splitting function p(qo,qi) = p(<7i,<7o) = p(qo)S(qo + qi — 2), we 
construct a related splitting function by breaking the go /<7i -symmetry: 

H(n „)-! 2 P(^,qi) (0<g <l<gi<2) 

p ^ qi >-\0 (0< 9l <l<g <2) . (46) 

The splitting function p(qo,qi) is almost identical to p{qQ,q\), their only difference being that the smaller splitting 
parameter qo is now always drawn for the branching into the left subinterval whereas the the larger splitting parameter 



1 The univariate version of this relation has already been discussed by Novikov in connection with the statistics of 
generalised multipliers, the so-called breakdown coefficients. 
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gi = 2 — go goes into the right subinterval. The one-point statistic is not affected by this modification of the splitting 
function, since it does not care whether the larger splitting parameter goes to the left or right subinterval. Hence the 
multifractal exponents are not changed either. With (^3|) we get 

f{v) = -v + — In ( / dqodqxpiqo, qi)(q% + <?i ) 



In 2 
1 



o 

2 



= -(? - 1) + ^2 ln (/ d( ZoP(<Zo)9o 
= rM ■ (47) 

We are hence forced to conclude that a multifractal analysis does not suffice completely to characterize the dynamics, 
i.e. the splitting function, of the underlying multiplicative cascade process. This statement holds even more strongly 
once we also give up energy conservation in the splitting function, with the consequence that properly defined backward 
moments show deviations from perfect multifractal scaling. In order to attain the full information, there is no way to 
get around spatial correlations in general, and the cumulant correlation densities Ck 1 ,....k ti in particular, to extract 
the branching moments ((lngo) n_m (hi<7i) m ) c and hence the branching generating function 

oo 1 

Q\>*M = £ -^((lng r(ln<7ir> c AW . (48) 
^ „ ni\n 2 \ 

ni ,^2—0 

Using Eq. ( p^ ) the branching generating function can then be uniquely inverted into the splitting function via a 
two-dimensional inverse Laplace transformation, 



/>oo 

/ dxdyp(2e- x ,2e- y )e 
Jo 



(X + l)x-(\ 1 + l)y _ Q[A ,Ai]-(A +Ai+2)ln2 



e yiA ,A 1 ]-(A +A 1 +2)ln2 _ 



As a consequence we view Q[Ao,Ai] as the natural and complete generalisation of the multifractal mass exponents 
t(v) for random multiplicative binary cascade processes. 

V. EXAMPLES 
A. Multiplicative binomial process: the p-model 



The p-model splitting function given in Eq. (|l|), when inserted into Eq. (22), determines the branching generating 
function Q[Ao, Ai]: 

Q[A , Ai] = In Q {cxp [A ln(l + p) + Aj ln(l - (3)} 

+ exp [A ln(l - /3) + Ax ln(l + /?)]}) . (50) 



We introduce the transformation 



Ao = ^(A+ + A_) , Ar = i(A+-A_) , (51) 



which leads to 



Q[Ao, Ai] = ln( \ /?2) A + + ln fcosh ( V ln 



I + (3 

2 ' ~ v V 2 V 1 -P 

Fig. 2(a) illustrates this branching generating function. 



(52) 
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3 

a) 



3 

b) 




FIG. 2. Branching generating function Q[\q, Ai], shown in the range —3 < Ao, Ai < 3, for (a) the p-model with f3 — 0.4, 
(b) the a-model with [3 — 0.4, (c) the energy-conserving SRST cascade model with /3 = 3.2, and (d) the energy-non-conserving 
SRST cascade model with /3 = 3.2. 

Now we calculate the derivatives of the branching generating function with respect to the conjugate variables Ao 
and Ai, which are needed for the cumulant densities Ck 1v ..,k„. For the first derivative we find: 



<9<3[A , Ai 



9A 



<9Q[A , Ai 



a=o 9A i 

ln(l - (?) 



1 



A=0 
1 



hi 



1 + 

1-/3 



tanh ( —A_ In ( — —— 
2 \l-{3 



A=0 



-Hi-P 2 ) 



(53) 



For an arbitrary derivative of Q[Xq, Ai] of order ni+ri2 > 1, we get the following relationship: 



i) 



ni+n2 



Q[Ao, Ai 







ni-\-n2 



Q[A ,Ai 



A=0 



(-1) 



A=0 

a n 1+ n 2 Q [AojAl 



d\ 



aA™ i+ " 2 



A=0 



(54) 
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Except for an alternating sign, the two-point branching moments ((lngo)™ 1 (higi)" 2 ) c are identical to the one-point 
branching moments ((lngo)™ 1+ " 2 )c; this only holds for the p-model and not for the other three models associated 
with the splitting functions (||)~(||). As a consequence of (54), we only need to calculate derivatives of Q[Xo, Ai] with 
respect to Aq. We use the intermediate step of (153) and write for n > 1: 



a"Q[Ao,A! 



1 



In 



A=0 



l + ( d^ 1 



(3 J \<9An -1 



tanh 



Since tanh a; is an odd function in x, i.e. 

tanh a; 



~3 

(<ry2m 



1 5_iL 7 

15 315^ 



-A In 



G2 



2835 



1 + /3 
1-/3 



(55) 



A o =0 



(2 2 " 1 - 1)5; 



rn — 1 



(2m)! 



2m 2m-l 



(56) 



with the Bernoulli numbers B2 m , all odd derivatives of the branching generating function Q[Ao,Ai] with respect to 
Aq vanish: 



U 



ii 
ii 

o 



II 

5 



u 




FIG. 3. Spatial cumulant densities for the p-model with j3 = 0.4 and J = 6: (a) second order Ck x ,k, 2 , (b) fourth order with 
three equal indices, Ck 1 =k 2 =k 3 ,k 4 , and (c) fourth order with two pairs of equal indices, Ck 1 =k 2 ,k 3 =k 4 - The indicated bin 
labels are related to the binary indices k = (fci • • • fcj) by k — 1 + y\ 
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d n Q[X ,Xi 



= 



(n = 3,5,7, 



(57) 



A=0 



The same conclusion could have also been obtained directly from ( |54| ) : for odd n > 3 we have d n Q[\o, Ax]/ <9Aq |a=o — 
— d n Q[Xa, Ai]/9A"|a=0j but from go/9i _ symmetry considerations of the splitting function we expect ((lngo)™)c = 
((lngi) n ) c ; hence, d n Q[Xo, Ai]/<9Aq |a=o = ((higo) n )c = 0. An immediate consequence of this result is that, for 
example, the spatial cumulant densities of third order vanish, CkiK^k^ = 0; see Eq. (|32|). 
For even derivatives of Q[Xq, Ai] with respect to Ao we find from (|55|) and d56|) that 



d n Q[X a , Ai 



oxs 



(2" - 1)B„ 



A=0 



hi 



1-/3 



(n = 2,4,6,...) 



(58) 



The results (|53j), (|54j), ( J57[ ) and ( |5§| ) can now be inserted into Eqs. (g0|)-(|3^) to determine the cumulant correlation 
densities Cm, ...,k„ = ((l ne Ki) ' ' ' (lne/c„))c- Fig. 3 illustrates the results for the two-point statistics: Ck x ,k 2 of second 
order, and Cn 1 =K 2 =K 3 ,K i and CK 1 =K 2 ,K 3 =K i of fourth order; note again that the third order Ck,i,k 2 ,k.3, vanishes. 
Figs. 4(a) and 4(b) show the second-order and fourth order cumulants as a function of the ultrametric distances 
respectively. 




FIG. 4. (a) Second-order cumulant density Ck 1 ,k 2 of the p- model (solid line) and a-model (dashed line) as a function of the 
ultrametric distance d,2 = dist(«i, K2)- Parameters values (3 — 0.4 and J = 6 were used, (b) Two projections of fourth-order 
cumulant densities for the p-model as a function of the ultrametric distance d^. 
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Since the p-model splitting function is energy-conserving and go -symmetric we can use relation ([l5]) to extract 
the multifractal exponents t(i>) from the branching generating function <5[Ao, Ai]. This leads to the well-known result 



It is illustrated in Fig. 5. Via 



t(u) = In 

v ' In 2 



d n r(y) 



dv r > 



v=0 



l-(3 



1 d n Q[X Q ,Xi} 



In 2 



(59) 



A=0 



= -^,i + ^((ln?o)") c 



(60) 



the derivatives of t{v) with respect to v are linked to the derivatives (p3[), ( |57| ) and (58) of the p- model branching 
generating function (p2|). 



15 \ 



10 



-5 



-10 




FIG. 5. Multifractal exponents t{v) of the p-model (dash-dotted line) with (3 = 0.4 and of the SRST cascade model (full 
line) with f3 = 3.2. 
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B. a- model 



The symmetric a-model is similar to the p-model except that it does not conserve energy in a cascade splitting. 
From the splitting function (||), its b.g.f. is found to be 



Q[A 0) Ai] - ±(\ + \i)Hl-f3 2 ) 



In | cosh 
■ In < cosh 



ft 



Ml> 



2 1 U 



■ft 



ft 



(61) 



clearly different from the b.g.f. (p2| ) for the p-model. Fig. 2(b) depicts the b.g.f. (^Tj). Note that for A = or Ai = 
the two expressions ([32]) and (|61[) become identical; consequently, the one-point branching moments of the a-model 
are identical to those of the p-model given in (p3j), (p7j) and (pjj): 



(62) 



/a n Q[A ,Ai] 




fd n Q{X ,X 1 ] 




V 5A» 


A— 0^ a — model 


I 9A« 





This is the reason why, in a multifractal approach, the two models look the same asymptotically. To see differences 
between them, one must consider two-point branching moments: for the a-model they all vanish, 



d n Q[X Q ,\i 



dX^dXl 



<(lng )" 1 (ln<Zi)"-" 1 ) c = ° . 



(63) 



A=0 



where 1 < m < n. We hence see that the two-point branching moments are sensitive to the violation of energy 
conservation in the splitting function. As a consequence of ( |62j ) and (63), the cumulant densities CK x ,...,K n of the a- 
and p-model now look slightly different. Fig. 4(a) compares Cn lt K, 2 of second order as a function of the ultrametric 
distance (p7|). 



C. SRST cascade 



The p-model is able to describe the multifractal aspects of the intermittent fluctuations occurring in the energy 
dissipation field in fully developed turbulence. However, due to its simplicity, the experimental multiplier distributions 
cannot be reproduced. In Ref. |l3[ , a modification of the p-model was proposed which accounts for the correct 
multiplier distributions; we call this modification the SRST cascade model. 

Insertion of its splitting function (^) into Eq. ( p2| ) yields an analytic expression for its b.g.f., 



Q[A ,Ai] =ln 



dzz^+^l-z)^- 1 



In 



r(2/3) 



(A + Ai)ln2 



(64) 



Since /3 — 3.2 and Ao ~ « Ai, so that Ao + ft > and Ai + ft > 0, the integral appearing in the first term of the 
right hand side can be identified with the beta function B(X + ft, A x + ft) = T(X + ft)T(Xi + ft)/T(X Q + Ai + 2/3). 
This leads to 



Q[A ,Ai] - (A + A 1 )ln2 + ln 



r(2/3) 



r(A + ft) r(Ai + ft) 



r(A + Ai + 2(3) T(ft) 



nft) 



(65) 



This result is illustrated in Fig. 2(c). 

Fr om (|6q ) the spatial cumulant densities Ck^,...,^ can be calculated in the straightforward manner presented in 
Sec. [II E[ In the lowest even orders, the results are very similar to those obtained for the p-model; the odd orders, 



however, are not equal to zero anymore. 

Making use of the relationship ( (45| ) the result (^5|) translates into the following multifractal exponents: 

, , , 1 . / T(2ft) T(u + ft)' 



(66) 
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Note that this expression, which is illustrated in Fig. 5, is only defined for v + [3 > 0; for i> + [3 < the integral in (J64 ) 
diverges. This is equivalent to the statement that the negative moments with v < —f3 of the splitting functions (p) 
and (0) do not exist. For the multifractal spectrum to exist over the full v range, however, moments of all orders, both 
positive and negative, must be finite. The absence of finite negative moments hence implies that it is not possible 
to construct full r(y) and f(a) curves for this specific splitting function. It is thus an example of a well-defined 
selfsimilar cascade process which cannot be described fully by the multifractal formalism. 



D. SRST cascade with no energy conservation 



Experimentally, the intermittent structures in the three-dimensional energy dissipation field of fully developed 
turbulence are observed on a one-dimensional cut. Although energy is conserved in three dimensions, this is probably 
not the case in one dimension. For the multiplicative branching models, this has the consequence that the splitting 
function p(qo,qi) cannot be expected to conserve energy. In this spirit, the expression (|]) represents an untested 
extrapolation of the SRST multiplier distribution (||), which has been deduced from one-dimensional data under the 
assumption of energy conservation |13| . 

Insertion of the splitting function (0) into Eq. yields the corresponding branching generating function: 



Q[A ,Ai] - (A + Ai)ln2 + ln 



/ T(2/3) r(A +/3) 
\T(Xo + 2P) T(/3) 



In 



T(2f3) r(Ai+/?) 
r(Ai+2/3) T(J3) 



(67) 



the explicit derivation is analogous to the one given in Sec. VC. The illustration of ( J67| ) is given in Fig. 2(d); it 
differs from the branching generating function (|6^) of the SRST cascade model with energy conservation. As in the 
p/a-model comparison, the one-point derivatives 



9"Q[A ,Ai 



<9AS 



A=0 



("EC) 



d n Q[X ,X 1 



dXJl 



A=0 



SRST 
(EC) 



(68) 



of the SRST cascade model with (EC) and without (no EC) energy conservation are identical, while for the two-point 
derivatives we find (1 < n\ < n) 



0"Q[A o ,Ai 



= 0^ 



A=0> 



SRST 
(noEC) 



d n Q[X ,X 1 



(69) 



A=0, 



SRST 
(EC) 



Consequently, the cumulant densities CKi,...,K n of the SRST-cascade model with and without energy conservation are 
different. 



VI. CONCLUSIONS 



With a clever change of variables from energy densities en to the singularity strengths an or lne^, we have derived 
an analytic expression for the multivariate generating function of binary multiplicative cascade models. The latter 
completely describes the n-point statistics, i.e. the spatial (cumulant) correlation densities of arbitrary order n. The 
key input has been a bivariate branching generating function, which is related to the underlying splitting function 
of the binary multiplicative cascade process via a two-dimensional Laplace transform. This branching generating 
function can be understood as a natural and, for selfsimilar binary cascade processes, complete generalisation of the 
multifractal mass exponents. While its properties completely fix the spatial correlation densities, the multifractal mass 
exponents do not. Various cascade models, relevant to fully developed turbulence, have been discussed to underpin 
this point. 

We have shown that, given that the experimentally measurable cumulants in (lne^ ,) are n-fold derivatives of 
the branching generating function, the latter can in principle be reconstructed from the former. With the help of 
Eqs. pq ) and (|4^), the b.g.f. can then be inverted into the splitting function via a two-dimensional inverse Laplace 
transform. In this way, the violation of energy conservation along one-dimensional cuts through the three-dimensional 
energy dissipation field can be inspected. 

Before this new approach is applied directly to (turbulence) data, however, a number of complications will have 
to be dealt with. In order to infer the branching generating function from the cumulant densities CK ± ,...,K n (to all 
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orders, in principle), a very effective representation of the latter has to be found; here, as in Refs. [ pL 0|Jl l| , a wavelet 
transformation might be useful to compress the information contained in the cumulant densities. 

Moreover, there is the problem of non-homogeneity: as a consequence of the hierarchical nature of the cascade 
evolution, the theoretical correlation functions are not invariant with respect to spatial translations, in contradiction 
to experimental measurements. It remains to be seen whether and in what way a scheme to restore homogeneity, as 
for example the one used in Jl6[ , influences or destroys the capability of inferring the branching generating function. 

Furthermore, the observed statistical dependence of multipliers [|l3 17 ll| has to be taken into account. Recent 
simulations |l^] indicate that this issue is closely linked to the restoration of homogeneity. This is in agreement 
with the conclusions reached by Nelkin and Stolovitzky |ll| by a different route, who argue that the experimentally 
proven dependence of multiplier distributions on the position of the subinterval implies that the multipliers are not 
statistically independent. Since any scheme to restore homogeneity will necessarily average out subinterval positions 
in some way, it will likely influence the multipliers' statistical dependence also. This remains to be explored in detail 

Finally, the binary structure of the selfsimilar cascade processes discussed in this paper may not be appropriate: 
assuming that the physical processes themselves are, indeed, selfsimilar, the best selfsimilar basis for a scaling analysis 
(such as a specific wavelet) should be selected by the data itself. 

Once these points are clarified, new information can hopefully be gleaned from the analysis of "fully developed 
turbulence data". Besides fully-developed turbulence, we envisage many and diverse applications of our analytic 
solution in other branches of physics. The case of QCD branching processes immediately comes to mind. For the 
latter, the a- and p-models have already been used in this context as simulation toy models [ p0|]2l| l. Implications in 
this and, for example, random multiplicative process calculations in large-scale structure formation in the universe 
p2[ remain to be explored. 
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